ENTROPY PRODUCTION IN THERMOSTATS II 
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Abstract. We show that an arbitrary Anosov Gaussian thermostat close to equi- 
hbrium has positive entropy poduction unless the external field E has a global 
potential. The configuration space is allowed to have any dimension and magnetic 
forces are also allowed. We also show the following non-perturbative result. Suppose 
a Gaussian thermostat satisfies 

K^(a) + i|£;,|2 <0 

for every 2-plane cr, where is the sectional curvature of the associated Weyl 
connection and E^r is the orthogonal projection of E onto a. Then the entropy 
production of any SRB measure is positive unless E has a global potential. A 
related non-perturbative result is also obtained for certain generalized thermostats 
on surfaces. 



1. Introduction 

In this paper we consider the dynamical system given by the motion of a particle of 
unit mass on a closed Riemannian n- manifold M subject to the action of an external 
field E. We also enforce as a constraint that the kinetic energy is a constant of motion, 
so the resulting equation is: 

where D denotes covariant derivative and 7 : R — > M is a curve in M. This equation 
defines a flow on the unit sphere bundle SM of M which reduces to the geodesic 
flow (free motion) when E = 0. The kinetic energy is held fixed by Gauss' principle 
of least constraint and thus the system defined by is referred to as Gaussian 
thermostat. Thermostats have become quite popular as models in nonequilibrium 
statistical mechanics [TTl HU [THJ EI] . Like geodesic flows, they are reversible, that 
is, the flip SM 3 {x,v) 1— >• {x, —v) G SM conjugates with (/)_t. 

Let G^; be the vector field in SM that generates 0. An easy calculation (cf. |,3 6. ) 
shows that the divergence divGe of G^; with respect to the canonical volume form 
O of SM is given by 

(2) divGij = -(n - 1)0, 

where 6 is the 1-form dual to E, i.e., Ox{v) = {E{x),v) and we regard 6 also as a 
function 6 : TM M. We see right away that (p does not preserve the Liouville 
measure (i.e. B) unless E = 0. But in principle, the flow may preserve other smooth 
measures. In fact, it is an exercise to check that preserves a smooth volume form 
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iff ^ is a coboundary, that is, iff there exists a smooth function u : SM M which 
solves the cohomological equation 



For example, suppose 6 is an exact 1-form, i.e., the external field E has a global 
potential U and write E = —VU. Then Ge{—U o tt) = 6, where vr : SM M is 
footpoint projection: it{x,v) = x. However, in general, one does not expect to have 
smooth solutions of 0- For example, if 6' is a closed, non-exact 1-form and every 
homology class in Hi{SM,Z) contains a closed orbit of 0, then there is no global 
solution to (jni)^ 

In the presence of hyperbolicity, there is a close relationship between Q and the 
entropy production of an SRB state p which we now describe. We will say that a 
0-invariant measure p is an SRB measure (or state) if p is ergodic and 



where hp{(j)) is the measure theoretic entropy of (p with respect to p. The entropy 
production of the state p is given by (cf. |29j) 



f^p{4>) = y, positive Lyapunov exponents = — negative Lyapunov exponents. 



Suppose now that (p is an Axiom A flow (we recall the definition in Section 3) and 
let p be an SRB state. We will see in Lemma (3.11 that if e^{p) = 0, then is in fact 
a transitive Anosov flow and (jH)) must hold. Conversely if (jH)) holds, then preserves 
a smooth measure and is a transitive Anosov flow. Hence p must be the unique 
invariant smooth measure and consequently (0)) holds which in turn implies e^{p) = 0. 

Thus, for Axiom A thermostats, e(^(p) = iff there exists a smooth solution of Q- 

In Section El we will explain why a transitive Anosov thermostat is always homo- 
logically full, i.e. every homology class in Hi{SM, Z) contains a closed orbit. Thus, if 
9 is closed, but not exact (e.g. electromotive forces), then 6^(0) > for any Axiom A 
thermostat. This was proved by M. Wojtkowski |35| Proposition 3.1] assuming that 
is an Anosov flow topologically conjugate to a geodesic flow, and by F. Bonetto, G. 
Gentile and V. Mastropietro PP for the case of a metric of constant negative curvature 
and 6 a small harmonic 1-form. 

The natural question now is: what happens for an arbitrary field E which does not 
necessarily have local potentials? In two degrees of freedom the problem was solved 
completely in [Hj: an Anosov Gaussian thermostat has zero entropy production iff E 
has a global potential. The aim of the present paper is to provide similar results for 
n degrees of freedom. 

^Note that n* : H^{M,R) H^{SM,R) is an isomorphism for M different from the 2-torus, 
thus 9 is exact iff Tr*9 is exact. 



(3) 



Ge{u) = 6. 




e^p^p) := — / divGEdp = — Lyapunov exponents. 
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We note that the assumption that (f) is uniformly hyperbohc is known in the htera- 
ture on nonequihbrium statistical mechanics as the chaotic hypothesis of G. Gallavotti 
and E.G.D. Cohen: for systems out of equilibrium, physically correct macroscopic 
results will be obtained by assuming that the microscopic dynamics is uniformly hy- 
perbohc. A system with e(/,(p) > is sometimes referred to as dissipative. Dissipative 
Gaussian thermostats provide a large class of examples to which one can apply the 
Fluctuation Theorem of Gallavotti and Cohen ^3 ^1 ^| (extended to Anosov flows 
by G. Gentile [15j) and this theorem is perhaps one of the main motivations for 
determining precisely which thermostats are dissipative. 

In our first result we will allow magnetic forces. This involves the addition of a 
Lor entz force F to the right hand side of (QJ). For each x G M, F^, : T^M T^M is an 
antisymmetric linear map such that the 2-form (F^(t;),w) is closed. We will indicate 
this thermostat by 0£;,f- Note that (/)o,f is a magnetic flow and hence it preserves the 
volume form 0. Suppose 0o,f is Anosov and E is an arbitrary external field. Then 
for e sufficiently small and s G (—£,£:), the flow (psE.Y is also a transitive Anosov 
flow. Moreover, the map [—e,e) 3 s y-^ e(s) := e0^^p(ps) is smooth jHUHOlEl]- It is 
immediate that e'(0) = and in Section 2 we will show that e"(0) > with equality 
iff E has a global potential. Thus we obtain: 

Theorem A. An Anosov Gaussian thermostat close to equilibrium has zero entropy 
production if and only if the external field E has a global potential. Magnetic forces 
are allowed at equilibrium. 

We now explain the non-perturbative results (which do not include magnetic forces). 
Given a 2-plane a C T^M, set: 

(5) k{a) := K{a) - div^E - \E\^ + 

where K{a) is the sectional curvature of the 2-plane a, Ea- is the orthogonal projection 
of E onto a and div^-i? := (V^-E, ^) + {Vr,E, rf) for any orthonormal basis {^, rj} of a. 
The expression 

K^{a) := K{a) - div^E - \E\' + \E„\\ 
is precisely the sectional curvature of the Weyl connection IH7] : 

V^Y = VxY + {X, E)Y+ {Y, E)X- {X, Y) E. 

Hence 

kia) = KUcT) + \\E^\^ 

Theorem B. Let (p be a Gaussian thermostat with k < and let p be an SRB 
measure. Then ep{(f)) = if and only if the external field E has a global potential. 

Like in [Hj this non-perturbative result will be established by using Pestov type 
identities as in IH IHI for geodesic flows. A closely related result about the cohomo- 
logical equation Ge{u) = where is an arbitrary 1-form is presented in Theorem 
Ol 
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We remark that in Theorem 5.1], M. Wojtkowski has shown that for n > 3, 
the condition k < imphes that (p is Anosov and for n = 2 it suffices to assume that 
Kyj < 0. (In general, < only ensures that the flow has a dominated splitting.) 

Our last non-perturbative result concerns a more general class of thermostats, but 
it will be only for n = 2. In principle, nothing impedes us from considering external 
fields acting on the particle which are also velocity dependent. The way to formalize 
this is to say that our external field is a semibasic vector field E{x, v), that is, a smooth 
map TM 3 {x, v) ^ E{x, v) e TM such that E{x, v) e T^M for all (x, v) G TM. As 
before the equation 

dt |7|2 
defines a fiow (j) on the unit sphere bundle SM. These generalized thermostats are 
reversible as long as E{x,v) = E{x, —v). 

Suppose now that M is a closed oriented surface. Set A(x, v) := v),iv), where 
i indicates rotation hj tt/2 according to the orientation of the surface. The evolution 
of the thermostat on SM can now be written as 

D'j 

(6) — = A(7,7)i7. 

If A does not depend on v, then is the magnetic fiow associated with the Lorentz 
force F^(f ) = \{x)iv. If A depends linearly on v, we obtain the Gaussian thermostat 

If we fix a Riemannian metric on M, its conformal class determines a complex 
structure. Given a positive integer k, let TCk denote the space of holomorphic sections 
of the k-th power of the canonical line bundle. By the Riemann-Roch theorem this 
space has complex dimension {2k — l){g — 1) for A; > 2 and complex dimension g for 
k = 1, where g is the genus of M. (For A; = 1 we get the holomorphic 1-forms and for 
k = 2 the holomorphic quadratic differentials.) Note that the elements in Tik can be 
regarded as functions on SM."^ 

Recall that vr : SM — > M is a principal S^-fibration and we let V be the infinitesimal 
generator of the action of S^ . If G denotes the vector field that generates the geodesic 
flow, the horizontal vector field H is given by the Lie bracket H = [V, G]. 

Theorem C. Let M be a closed oriented surface and consider an Anosov generalized 
thermostat (0) determined by X = 3?(g), where q G Hk- Suppose 

K - H{X) + X^[{k + 1)V(2A; + 1)] < 0, 

where K is the Gaussian curvature of M. Then has zero entropy production if and 
only if X = 0. 

When K = —1, k = 1, and A is sufficiently small, the theorem is proved in [I] using 
the same perturbative methods as we will use for the proof of Theorem A. Note that 
for k odd the flow is reversible, so Theorem C provides a large class of new examples 
to which the Fluctuation Theorem of Gallavotti and Cohen applies. 

^Sections of the fc-th power of the canonical hne bundle can be regarded as functions on SM 
which transform according to the rule f{x,e^''^v) — e'^'"^ f {x , v) . 
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2. Derivatives of entropy production 

2.1. The variance. Let be a transitive Anosov flow on a closed manifold X. We 
will assume that is weak-mixing, i.e., the equation F o 0j = e*'**-F, a > 0, has no 
continuous solutions. 

Let yU be a Gibbs state associated with some Holder continuous potential. Given a 
Holder continuous function F : X — > R, the variance of F with respect to /i is deflned 
as: ^ 

Var^(F):=^lim (Fo0,-F)rft^ rfp, 

where 

Jv 

This limit exists and it appears in the central limit theorem for hyperbolic flows [27 \ . 
There are other equivalent ways of expressing the variance. Let 



PF{t) := [ {Foct>fF-F^)dii 

JX 



be the auto-correlation function of F . Then the variance can also be expressed as (cf. 
Section 4]): 



pp{t)dt = 2 pF{t)dt. 

-oo Jo 

In fact the Fourier transform of pi? 

pUw):= / e'^'pF{t)dt 



deflned as a distribution, has a meromorphic extension to a strip < e with no 

pole at w = j2Hll2ni- The value at w = is precisely Var^(F)/2. 

2.2. Proof of Theorem A. We flrst recall the setting described in the introduction. 
Consider a closed Riemannian manifold and F a Lorentz force. Suppose 0o,f is Anosov 
and E is an arbitrary external fleld. Then, by structural stability, for e sufficiently 
small and s G (—£,£), the flow (j)sE,v is also a transitive (and weak-mixing) Anosov 
flow. 

Consider the map 

(-£,£) 3 e(s) := e^^^^{ps). 
It follows from the results of G. Contreras [3] or Ruelle jHOJ E21 that this map is 
smooth. Indeed since divGs£;,F = — ~ we have 



e(s) = (n — l)s / 9 dp, 

J SM 



and the results in jHUl E2] assert that s ^ J^^j 9 dps is smooth. Thus 



e'(0) = (n- 1) / 9 dp, 

JsM 
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e"(0) 



9dps. 



s=0 -JSM 

Here p := po is the Liouville measure of SM. We see right away that e'(0) = since 

SM 

because 9x{v) = —9,j.{—v). The derivative 

d 
ds 



Odps 



s=0 JSM 

can be computed from the results in [SOIIS^]- Given a smooth function F : SM 
the derivative 

■ Fdp, 



s=0 JSM 

is the hmit as u ^ with Q(w) > of 



SM 



Awt 



d{F o (y (x, v)) dp^x, v) 



divF(a;, v) F{(j)t{x, v)) dp{x, v) 



SM 



where Y is such that Gse,f = Go,f + sY. Since div (Gs£;,f — Go,f) 
see that 



-s {n — 1)6 we 



e"(0) = 2{n - i f hm 



Awt 



6^{v) 6{4>t{x, v)) dp{x, v). 



SM 



As pointed out before 



Awt 



9^{v) 9{(f)t{x,v)) dp{x,v) 



SM 



extends to a holomorphic function near w = 0. We can now identify 



2 hm 

w^O 



Awt 



Ox{v) 9{(f)t{x,v)) dp{x,v) 



SM 



with the variance of the function 6 with respect to the LiouviUe measure p. Thus 

e"{0) = {n-lfYa.T,{e). 

The variance has the wonderful property that Var^(F) > with equality iff F is a 
coboundary ([23 Section 4]). Thus we have shown that e'(0) = and e"(0) > with 
equality iff there exists a smooth solution u to the cohomological equation 

(7) Go,F(n) = 6. 

But the results in ^ Theorem B] give a complete understanding of the cohomological 
equation for Anosov magnetic flows. Indeed, there is a solution of ((Tj) iff 6' is an exact 
form. For geodesic flows (i.e. F = 0) this result is proved in |8^. 

Thus, unless E has a global potential, e"(0) > and therefore e(s) is strictly 
positive for s ^ near zero. This shows Theorem A. 
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2.3. Some explicit calculations of Va,T^{6). Suppose M is a compact locally sym- 
metric space of negative curvature and suppose 6' is a harmonic 1-form. If F = 0, 
then Var^(^) has been calculated by A. Katsuda and T. Sunada in PTE Proposition 
1.3]. They show that 

where h is the topological entropy of the geodesic flow of M. For n = 2 and h = 1 
{K = —1) we obtain 

With an appropriate normalization for the L^-norm of 6 we recover precisely the 
calculation performed in [T} Page 687] to compute e"(0). 

It is interesting to see what happens for n = 2 and K = —1 if one adds a uniform 
magnetic field. Suppose we take F^^v) = iv and A G [0, oo). It is well known that for 
< A < 1, the fiow 0o,AF is Anosov and for A = 1 we obtain the horocycle fiow. Let 6 
be a harmonic 1-form and for A G [0, 1) let us try to compute VarA,^(^), the variance of 
9 with respect the fiow 0o,af and the Liouville measure /i. A direct calculation along 
the lines in is possible, but we will take a different, more economical approach that 
exploits the good properties of the variance. It is also well known (see for example 
|22j) that the flow 0o,af is conjugate to the geodesic flow 0o,O; up to a constant time 
scaling by a/I — A^. Let f = f\ : SM — > SM be this conjugacy and note that it is 
immediate to check that /o is the identity, so / is isotopic to the identity. Since / is 
a conjugacy: 

where = Vl — A^ Go,o- Observe that Ox{v) = 'k*9{Go^xf){x,v) and therefore 
9 o f{x,v) = 7r*9{Go,xFKf{x,v)) = n*9{df^,,,)Xx) 
= r7r*9{Xx){x,v). 

Hence 

VaixAO) = Varx„;.(^ o /) = Varx„;.(/*vr*^(XA)). 

Since 9 is closed, f*n*9 is a closed 1-form in SM. Observe that Varx;^,^(/*vr*6'(XA)) 
only depends on the cohomology class [f*7r*9] since the variance vanishes on cobound- 
aries. We noted before that / is isotopic to the identity, thus 

Varx„^(/*7r*^^(XA)) = Varx„^(7r*^^(XA)) = YaTj^„^{VT^ 9) 

= (l-A2)Varx„,W. 

From the deflnition of the variance it follows right away that if F : SM — > M is any 
function then 

Varx„;.(F) = Vl^Va.T^{F) 

which yields 

Varx,,.(/*^*^^(XA)) = (1 - A^)^/' Var^(0). 
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Summarizing 

Var,,^(0) = (l-AY/'Var^(^). 

Thus we have obtained the following formula for the second derivative of entropy 
production in the presence of a uniform magnetic field with intensity A: 

A completely analogous formula can be obtained for compact quotients of complex 
hyperbolic space with magnetic field given by the Kahler 2-form. 

3. NON-PERTURBATIVE RESULTS 

All the results in the this section will be based on studying the cohomological 
equation using Pestov type identities. The results on the cohomological equation are 
all collected in Section^ 

3.1. Axiom A thermostats. A closed 0-invariant set A is said to be hyperbolic if 
T{SM) restricted to A splits as Ta{SM) = RGe ® E"" ® in such a way that there 
are constants C > and < p < 1 < r/ such that for alH > we have 

< Cr/"* and < Cp*. 

The flow is Axiom A if the nonwandering set f2 is hyperbolic and the closed orbits 
are dense in VL. Recall that by the Smale spectral decomposition, is a finite union 
of disjoint basic hyperbolic sets. A hyperbolic basic set is a hyperbolic set such that: 

• the periodic orbits of 0|a are dense in A; 

• 0|a is transitive; 

• there is an open set [/ D A such that f\tm4'tiU) = A. 

The flow is Anosov if SM is a hyperbolic set. If is Anosov, it is also Axiom A 
(but not conversely, of course). Recall that there are examples of Anosov flows for 
which Q is not the whole space j^]. 

Lemma 3.1. Let be an Axiom A thermostat and p an SRB state. If ep{(f)) = 0, 
then (f) is a transitive Anosov flow and there exists a smooth solution u of Ge{u) = 9. 

Proof. Let A be the basic hyperbolic set on which p is supported. Since p is an SRB 
measure, [21 Theorem 5.6] implies that A is an attractor, that is, there exists an open 
set D A such that A = r\t>Q(t)t{U) . 

As pointed out in the introduction, if 6^(0) = 0, then 

hp{4>) = positive Lyapunov exponents = — negative Lyapunov exponents. 

Thus p is an SRB measure for both 0^ and 0_t. Consequently, A is an attractor for 
both 0i and 0_f. This forces A to be open and since it is closed, A = SM and is a 
transitive Anosov flow. 

Let J* and J"^ be the stable and unstable Jacobians of 0. If p is an SRB measure for 
both 0t and 0_t then the theory of Gibbs states for transitive Anosov flows (cf. |19| 
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Proposition 20.3.10]) implies that — ^j^^^log J^" and -^l^^^^og are cohomologous 
(and the coboundary is the derivative along the flow of a Holder continuous function). 
It follows that preserves an absolutely continuous invariant measure with positive 
continuous density (and this measure would have to be p). An application of the 
smooth Livsic theorem fIT\ Corollary 2.1] shows that preserves an absolutely con- 
tinuous invariant measure with positive continuous density if and only if preserves 
a smooth volume form. But if preserves a smooth volume form, then there is a 
smooth solution m of G£;(m) = 6^ as desired. 

□ 

3.2. Proof of Theorem B. We are required to prove that if ep(0) = 0, then E has 
a global potential. By [37, Theorem 5.1] the condition A; < implies that is Anosov 
and hence Axiom A. By Lemma l3.ll there is a smooth solution u of Ge{u) = 9. 
Theorem 14.41 implies that 9 is exact as desired. 

3.3. Proof of Theorem C. We will need some preliminaries which can all be found 
in [TT]. Let L'^{SM) be the space of square integrable functions with respect to the 
Liouville measure of SM. The space L'^{SM) decomposes into an orthogonal direct 
sum of subspaces ^ Hn, n E 1^, such that on —iV isn times the identity operator. 
Consider the following first order differential operator: 

r/_ := {G + zH)/2. 

The operator ?7_ extends to a densely defined operator from Hn to Hn-i for all n. 
If we let C^{SM) = Hn n C^{SM), then r/_ : ^ C^_^ is a first order elliptic 
differential operator. The kernel of the elliptic operator ■)]_ in C^{SM) is a finite 
dimensional vector space which can be identified with Tik. (For all these properties 
see tlli.) 

Now take q G Hk and let A := 3?(g). Then p := V{\) = ^{ikq). Since ri^q = 0, we 
see right away that G{p) + HV{p)/k = and hence Theorem 14.61 implies that V{X) 
is a coboundary iff A = 0. But ^(A) is the divergence of the generalized thermostat 
with respect to 6 (cf. Lemma 3.2]), so the entropy production vanishes iff A = 
as desired. 

4. PESTOV IDENTITY AND COHOMOLOGICAL EQUATION FOR THERMOSTATS 

4.1. Semibasic tensor fields. Let vr : TM \ {0} — M be the natural projection, 
and let := 7r*rJM denote the bundle of semibasic tensors of degree (r, s), where 
rJM is the bundle of tensors of degree (r, s) over M. Sections of the bundles 
are called semibasic tensor fields and the space of all smooth sections is denoted by 
C°°(/3JM). For such a field T, the coordinate representation 

holds in the domain of a standard local coordinate system {x\y^) on TM \ {0} 
associated with a local coordinate system (x*) in M. Under a change of a local 
coordinate system, the components of a semibasic tensor field are transformed by the 
same formula as those of an ordinary tensor field on M. 
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Every "ordinary" tensor field on M defines a semibasic tensor field by the rule 
T 1-^ T o TT, so that the space of tensor fields on M can be treated as embedded in 
the space of semibasic tensor fields. 

For a semibasic tensor field (Tj^^ ' j^)(x, y), the horizontal derivative is defined by 

d ■ ■ d ■ 

3l...js\k ~ Q^k ji---js kqiJ Qyp ji...js 



_|_ \ ^ ■pimrpil---im-ipim + l---ir \ ^ "pP rpii . . .i 

kp-'-ji-.-js 2-^ kjm 

and the vertical derivative by 



m-lPJm + l---Js' 

m=l m=l 



d 

The operators 

V| : C^iPlM) C^{(3l^^M) and V. : C^{(31M) C^{(3l^^M) 
are defined as 

(V|T)-:-, = V|.r-:t := t-;-,, and (v.T)-:-, = v..i-:.j: = r-;-.,. 
/i jj 

In 1231122], the operators V| and V. were denoted by V and V respectively. 

4.2. Thermostats and the modified horizontal derivative. Let {M,g) be a 
closed connected Riemannian manifold. Given a vector field E on M, define Y G 
C^iPlM) by 



We have seen that the equation 



\y\ 



^ = >^(7,7)(7) 



defines the Gaussian thermostat on SM. 

Since the flow on SM defined by (jHI) depends only on the restriction of Y to SM, 
we may redefine Y to be 

without changing the flow on SM, so that from now on 

Ylix,y) = y,E\x)-E,{x)y\ 
Given T = (Tj^\" j^) G C°°(/3JM), we define the modified horizontal derivative as: 

r-pil...ir _ rpil...ir i rpii...ir yj 

so that V: : C~(/3JM) ^ C°°(/3J+iM). 

For convenience, we also define V', V', and V' as 
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We also set 

..jsik- 



-^ijiii...ir _ krpii...ir 



In particular, if u G C°°{TM \ {0}), then 

Note that X restricted to SM coincides with Gg. 
It easy to see that if 7 satisfies (jHI), then 

(Xr)(7,7) = ^(T(7,7)). 
Straightforward caluculations give: 

(9) Y;^, = y,E% - E,,,y^ 



(10) Yj., = g,,E^-E,5l 



(11) y'k = yi = VkE' - E,y\ 

where (y*) stands for the semibasic vector field (x, ?/) t-^ (y*). 

For V = {V') G C^{(3^M), we set 

h V m 

div V := V\l, div V := V.l div V := V.}. 

We recall the Gauss-Ostrogradskii formulas for the horizontal and vertical diver- 
gences inSI (see also Section 4.2], which deals with the case of Finsler metrics). If 
V{x, y) is a smooth semibasic vector field positively homogeneous of degree A in y, 
then 

(12) [ divVdfi = 0, 

J SM 

(13) / divVdii= I {X + n-l){V,y)dfi, 

J SM J SM 

where is the Liouville measure on SM . Whence we also have 

(14) f d\YVd^={\ + n) [ {E,y){V,y)dfi-{\ + l) [ {V,E)dfi, 

JsM J SM JsM 



because 



dTv F = = + V^,{ykE' - Euy') = div + div ((\/, y)E) - (A + 1) {V, E). 
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4.3. Pestov identity. Given a function u : SM —>■ M, we will also denote by u its 
extension to a positively homogeneous function on TM \ {0} (hoping that this will 
not yield any confusion). 

We first recall commutation formulas for horizontal and vertical derivatives [SH! 
(see also [7| Lemma 4.1], which deals with the case of Finsler metrics). If m e 
C°°(TM\ {0}), then 

(15) u.ik - u.k-i = 0, 

(16) uii.k - u.k\i = 0, 

(17) u\i\k - u\k\i = RlkU.i, 

where R},^ = y^R)ik and R is the Riemann curvature tensor. 
The next lemma is an analog of [H Lemma 4.5]. 

Lemma 4.1. If u e C°°(TM\ {0}), then 

(18) ua-k - u.kd = {gikE' - Ei6l)u.i, 

(19) U-Mk - U;k:l = RlkU-i, 



with 

Proof. We have 

whereas 
Thus, 



^Ik — ^kl + 0^k\l ~ ^|fc) + O^k^l-j ~ ^l^k-j)- 

Ud-k = {u\i + Ylu.i).k = u\i.k + Yl^u.i + Y^u.i.k, 
u.kd = u.k\i + Ylu.k.i- 



u.. 



ki-k - u.kd = {u\i.k - u.k\i) + Yi\u.i + Yl{u.i.k - u.k.i). 

Using (HSl), (USD, and (HUD, we come to (fTH|l . 
Further, 

u-d-.k = u,i\k + Y^u.d.j = (m|; + Y{'u.j)ik + Y^{uii + Y'u.s).j 
Therefore, 



I'.i.k - u,kd = {u\i\k - u\k\i) + (y^j^, - Yj!^^)u.j + Yi\u.j\k - u\k.j) 



u. 



k jU.s-j- 



In view of renaming indices and regrouping, we come to (jl9|) . □ 
The next lemma shows a Pestov type identity for thermostats. 



ENTROPY PRODUCTION IN THERMOSTATS 11 13 

Lemma 4.2. If u E C°°(TM\{0}) is homogeneous of degreeO iny, then the following 
holds on SM : 

m V 

(20) 2{V-u, V(Xm)) = I V^Mp + X((V M, V-u)) - div ((Xm)V m) + div ((Xm)V^m) 

- (Rj,(Vm), Vm) - (E,y)(VM, V^m) - (n- 1)(Xm)(E, Vm). 
Proof. With the above notations, we can write 

Xm = 

Therefore, 

(21) 2(V.(Xm), V.u) - div ((Xm)V^m) = 2g'^(Xu).iU.,j - {{'Xu)g'^u.,j).^ 

= g'^{Xu).^u.,j - {'Xu)g'^u.,j.i = 1-11. 
We rewrite the first term on the right-hand side of fIT\i as follows: 

I = 9'^{y''u:k);U;j = g'^{u;i + y^U;k-i)u;j 
= g'^u,iU.,j + g'^y^[u.i.,k + {u.,k.i - u.i.,k)]u.,j 

= iV^Mp + y^ig'^u.iu.,.).,k - y^g'^u.iu.,j.,k + g'^y''{gkiE"' - Ek5J')u.^u.,y 

Next 

y\g''u.,u..^)..k = M{^n,Vu)), 

y^g'^u.iu.,j.,k = y''g'^u.i[u.,k;j + {u.,j.,k - u,k:j)] 

= g'^u.Xy^u.,^):j - g'^u.iy^jU.,k + y'^g'^u.iKflu.m 

= (Vu, V^(Xu)) + (Xu)(E, Vu) + {ilyiVu), Vu) 

because 

g'H.iy^jU.,k = g'^u.iY^u,k = y'u.iE^u,k - E'u.iy^u,k = -{Xu){E, Vm), 

and 

g'^y\gkiE"^ - Ek6^)u.mU:j = y^E^'u.mU-.j - g'^y^EkU.iU.,j 

= {Xu){E,Vu) - {E,y){Vu,V-u). 

Thus, 

(22) / = \V-u\^ + X((Vm, V-u)) - (Vm, V^(Xm)) 

- {Ry{Vu),Vu) - {E,y){Vu,V-u). 
We rewrite the second term on the right-hand side of (PT|) as 

// = {Xu)g'^U.j.i = (Xu)g'^[u.i:j + [U-j.i - U.i:j)] 

= [iXu)g'^u.,]., - iXu)..,g'^u., + iXu)g'\g,,E"' - E,5T)u.^. 

Note that 

m 

[{y.u)g'^u.i\..j = div((XM)VM), 
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that 

{y.u).J^u., = {Vu,v\y.u)), 

and that 

{y.u)g'^{g,iE^ - E,6r)u.m = {n - l){-Xu){E,Vu) . 

Thus, 

m 

(23) 11 = div {(Xu)Vu) - (Vm, V^(Xm)) + (n - l){Xu){E,Vu). 
Inserting (j221)-(|2Sl) in dH)), we come to □ 
Note that for the curvature term in pO|) we have, putting Z = Vu: 

(24) {Ky{Z), Z) = (R,(Z), Z) - {VzE, Z) - (E, Z)\ 
Indeed, 

(R,(z), z) = [R,, + (r,^|, - + (F,^■l^^ - F/r,^)] y^z^z, 

= (R,(Z), Z) + (i/fcEj - - t/^E;, + Ei^^y')y'Z^Zi 

+ [(yfcE^' - E,,y^){gi,E' - E^d^) - {y^E^ - Eiy^){gkjE' - EJ^]y'Z^Z, 

= (Ry(Z), Z) — E'^f^Z^Zi — EkEiZ'^Z\ 

where we used the fact that {Z, y) = by homogeneity. 
One more useful identity is: 

(25) X(Vm) = V (Xu) - V'-u - {E, Vu)y + {E, y)Vu. 
Indeed, 

X(m-^) = y''{g'^u.j).,k = i/V^(M:fc.j - {u.,k-j - u.j.,k)) 

= g'Ky^n..k)., - g'^u., - g'^y\gk,E^ - Ek5f)u.^, 

and since 

g'^y\gk,E^ - EJf)u.^ = {E,Vu)y' - {E,y)u\ 

we have ((23) • 

4.4. Cohomological equation. Suppose that the cohomological equation 

(26) Geu = ^ 

holds with a smooth function u on SM and a smooth 1-form ^ on M. Denoting the 
homogeneous extension of u to TM \ {0} by -u as before, we get 

X.u{x,y) = {F{x),y), 

where F is the vector field dual to 'd with respect to the Riemannian metric. 
Integrating against the Liouville measure dfi and using (|T!?jl . (fT^ yields 

2 [ {Vu,V{Xu))dfi= [ \\V-u\'^ + X{{Vu,V-u)) + n(Xuf 

J SM J SM 

- (Rj,(Vm), Vm) - {E,y){Vu,V-u) - {n - 2){Xu){E,Vu)\ d^. 
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Note that by (0) 

/ X{{Vu,V'-u))diJ={n-l) {E,y){Vu,Vu)dn. 

JsM JsM 

Therefore, 

(27) 2/" {V-u,W{Xu))dfi= [ \\V-u\^ + n{Xu)^ - {Ry{Wu),Wu) 

JsM J SM 

+ (n - 2) [{E, y) {Vu, Vu) - (Xm) {E, Vu)] } rf/x. 

Using (P^. we have 

(E,y)(Vn,V^M)-(Xn)(i?,Vn) 
= -(E,y)(VM,X(VM)) + (E,y)(Vn, V(Xn)) + (E, y)2(V n, V^/) - (Xn)(E, Vn) 

= -(E, y) {Vu, X{Vu)) + {E, yy\Vu\' + div {u{E, y)V{Xu)} - u{E, V(Xn)) 

- u{E,y) dw {V(Xu)} -div {u(Xu)E} + u{E,V{Xu)) 
= -{E,y){Vu,X{Vu)) + {E,yf\Vu\^ 

+ div {u{E,y)V{Xu) -u{Xu)E} -u{E,y) div {V {Xu)}. 
Plugging this in (jTfjl and again using (fT!^ . we derive: 

(28) 2/ {V-u,V iXu))dfx= [ \\V-u\'^ + n{XuY - {'Ry{Vu),Vu) 

JSM J SM ^ 

in - 2){E,y){Vu,X{Vu)) + (n - 2){E,y)^\Vu\^^ dfi, 



where we used the equahty div [V (X'u)] = divF = 0. 
Since 



/ {\V{Xu)\'~n{Xur]df,= [ {\E\'-n{E,yf]d^i = 0, 

J SM JSM ^ 



' SM ^ ^ JSM 

we can rewrite (I^Hl) as follows, with Z = Vu: 



2/ {Vu,F)dfi= [ {\V-u\' + \F\'-in-2){E,y){X{Z),Z) 

JSM JSM 

-{Ry{Z),Z) + {n-2){E,y)^\Z\^]dfi, 



or 



(29) / ||F - Vu\' -in- 2){E,y){XiZ), Z) 

JSM 

- (R,(Z), Z) + (n - 2){E, yf\Z\^] d^i = 0. 
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From we obtain 

{-K{Z),Z) = {F - Vu,Z) + {E,y)\Z\\ 
Fixing any real parameter a, we now rewrite (j^^ as follows: 

/ { |F - Vu\^ - 2a{E, y){F - V'u, Z) - {n - 2 - 2a){E, y){X{Z), Z) 

JSM ^ 

- (R,(Z), Z) + {n-2- 2a){E, yf\Z\^] dfi = 0, 



or 



(30) 



f J\F-Vu- a{E, y)Z\' - (n - 2 - 2a){E, y){X{Z), Z) 

JsM 

- {tly{Z), Z) + (n - 2 - 2a - a^){E,yf\Z\'^'^ dfx = 0. 

Notice that 

2{E,y){X{Z),Z) = X{{E,y)\Z\') -Xi{E,y))\Z\^ 
A direct calculation gives 

X{{E,y)) = {\/yE,y) + \E\' - {E,y)', 

whence 

2{E,y){X{Z),Z) = X {{E,y)\Z\') - {{VyE,y) + \E\' - {E,y)')\Z\^ 
and therefore 

2/ {E,y){XiZ),Z)df, 

JsM 

= [ X{{E,y)\Z\') df,- [ {{VyE,y) + \E\'-{E,y)')\Z\'d^ 

JSM JSM 

= [ \n{E,y)'\Z\' - {VyE,y) -\E\']df^. 

JSM ^ ^ 

Plugging this in (IHUj). we get 



/ {\F- Vu - a{E, y)Z\^ + ^ \ [(V,i?, y) + \E\^] 

JSM ^ ^ 



(R,(Z),Z)- 



n-2-2a\ fn-2 



2 



{E,y)^\Z\^]d^^ = Q. 



Changing " ^° 1— >• a and using (|^. we deduce: 



(31) / I |F - V^^i - (n/2 - 1 - «) (E, + « [(V,E, y) + 

- (R,(Z), Z) + (Vzi?, Z) + (E, - [a^ + (n/2 - 1)^] (E, (i/x = 0. 
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So, if 

i^K,) - {V^E, - a{V,E, r,) - - (E, if + [a^ + (n/2 - 1)^] {E, r,f < 

for every x G M and every pair of orthogonal unit vectors ^iTj & T^M, then Z = 0. 
This means that m is a hft to SM of a function on M, u{x,y) = <f{x), and the 
cohomological equation imphes: 9 = dip. Choosing a = 1 and putting 

ki{a) = K{a) - diY^E - \E\^ + [l + (n/2 - 1)^] \E^\^, 
n .\2,„ ,2 



we arrive at the following: 

Theorem 4.3. Suppose ki <0 and Ge{u) = d. Then -d is exact. 

It is interesting to notice that ioi n = 2, ki = K — div E = K^, and that for n = 3, 
ki equals k of 

4.5. Using the invariant measure. The measure f/i is invariant if G£;(log/) = 
{n — 1)6. We let u = log/, so that 

(32) Ge{u) = {n-l)e 

and u = e"/i is a flow invariant measure. 

Let V{x, y) be a smooth semibasic vector field positively homogeneous of degree A 
in y. Since 

e"divl^ = e"Ki = div (e"l^) - {V, VM)e", 

we have by (fTSj) 



(33) 



I divVdu= [ \{X + n-l){V,y) - {V,Vu)]du, 

J SM J SM ^ ^ 



and, since 

m m 

e'^diYV = e^V^l = div (e"y) - (V, V-u)e'', 

we have by (jl4p 

(34) [ (SvVdu= [ \{X + n){E,y){V,y) - {X + l){V,E)du- {V,Vu)]du. 

J SM J SM ^ 

Let us integrate (^01) against dv. Using the flow invariance of u together with 
and dSll) yields: 

(35) 2 f {V-u,V(Xu))du= [ \\V-u\'^ + n{Xuf - {n - 2)(Xu){E,Vu) 

JsM JsM ^ 

- {E, y) (Vu, V-u) - (Rj,(V m), Vm) } dv. 
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We have 

n I {y.ufdv = n{n-lf! {E,y)'^e''dfx={n-lfn[ {E , y) {e"" {E , y)) dfx 

JsM JsM JsM 

I' 

\2 / j;,,/'„w/77i „.\TP\J,. /„ T\2 / //771 „,\/V7-„, ZT'X I I ZPISn 



= in-iy / divie^{E,y)E)dfi=in-iy / {{E,y){Vu, E) + \E\') du. 

JSM JSM 

Plugging this, ipij) and (jH^ in we have, with Z = V m: 

2(n-l) /" {V'-u,E)dv= [ \\V-u\^ + {n-l)e{E,Z) + {n-lf\E\^ 

JSM JSM 

- e{Z, V-u) - (RyiZ), Z) + {VzE, Z) + {E, Zf^ dv = 0, 



or 



(36) / {\Vu-{n-l)E-{l/2)eZ\'-\e'\Z\' 

JSM ^ 

- (Rj,(^), Z) + {VzE, Z) + (E, rfz/ = 0. 



So, if 



for every a; G M and every pair of orthogonal unit vectors C,,"!] ^ T^M, then Z = 0. 
Passing by we note that this condition also implies that is Anosov by jSEl Theorem 
4.1]. 

Using (f25|l . we have 



(37) X{Z) = {n-l)E -V-u - {E,Z)y + eZ. 

Then we can rewrite (I36p as 

/■ {\MZ) + {E,Z)y - {i/2)ez\' -\e'\z\' 

J S Kd 

- (R,(Z), Z) + (Vz^, ^) + {E, Zf^ du = 



or 



/ { |X(Z) + (E, Z)y + (1/2)^Z|2 - 2^^(X(Z), Z) - 

- (R,(Z), Z) + {VzE, Z) + dzy = 

or, using 

2^^(x(z), z) = x(^^|zp) - \z\^^e = x{e\z\') - {{VyE, y) + |Ep - e^)\z\\ 
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as 

(38) / \\X{Z) + {E,Z)y + il/2)eZ\' 

JSM 

- {RyiZ), Z) + (Vz^, Z) + (V,E, y)\Z\^ + \Ef\Z\^ - {E, Zf - ^e^\Z\^^ du = 0. 
Recall that 

k{a) = K{a) - div^E - \E\^ + ^I^^T- 

Hence if k{a) < for every x and every two-plane a G T^M, (j38p implies Z = and 
hence we obtain: 

Theorem 4.4. Suppose k < and Ge{u) = d. Then 9 is exact. 

To complete these results we now show that if the thermostat is transitive and the 
inequality k{a) < holds for all a, then Z = 0. Indeed, in this case 

K{a) - div^E - \E\'^ + ^{E,r])^ < 

unless (E,^) = 0, where {r],^} is an orthonormal basis of cr. 
Then §^ yields 

{E,Z) = 

and 

XZ+{E,Z)y + {l/2)eZ = 0. 

Then 

XZ = --0Z, 

yielding 

X(|ZH = -9\Z\^. 
Assuming the set {Z 7^ 0} to be nonempty, we obtain on this set 

X(log = -0 = -{Xu)/{n - 1). 
Consider this equation on a dense orbit. Then 

on every connected component of the intersection of this orbit with the set {Z 7^ 0}, 
with some nonzero constant C depending on the component. Such a component is 
obviously open in this orbit. At the same time, it is closed as the right hand side of 
the above equality is separated from zero. This means that the whole orbit is in the 
set {Z 7^ 0} and so the above holds on this orbit with the same nonzero constant, 
which means that |2^| is separated from zero on a dense orbit, and consequently it is 
nonzero everywhere. We now show that this is not possible. 

Recall that V u := (u'*) where m'* := g'^^u.j and u.j := Fix Xq E M and consider 
the restriction -u of m to SxM. Since SxM is compact there is yo G S^M for which 
dyf^u = 0. Since u is homogeneous of degree zero we must have V m(xo, yo) = 0. 

Thus Z = Vu = everywhere in SM. 

Summarizing, we have proved: 
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Theorem 4.5. Let (p be a transitive Gaussian thermostat with k < 0. Then pre- 
serves a smooth volume form if and only if E has a global potential. 

4.6. Cohomological equation for generalized thermostats on surfaces. Con- 
sider the thermostat determined by an arbitrary function A G C°°{SM) and let Ga 
be its infinitesimal generator. 

Theorem 4.6. Let p G C°°{SM) be such that G{p) + HV{p) /k = for some positive 
integer k, and suppose 

K - H{\) + \''[{k + 1)V(2A; + 1)] < 0. 
Then there exists u G C°°(5'M) such that Gx{u) = p if and only if p = 0. 

Proof. Note that Ga = G + XV. We will use the following Pestov type integral 
identity proved in [HI Equation (13)]. Given u G C^{SM) we have: 



(39) 2 / HuVGxudii= / {Gxuydfi+ / {HuY d^ 

'SM JsM JsM 



[ {K - H{\) + \^){Vufd^i. 

J SM 



Using that G(p) + HV{p)/k = and that H and G preserve the Liouville measure 
we obtain: 



HuV{p)dfi = — / uHV{p)d^ = k / uG{p)d^ = —k I G{u)pd^. 

SM JSM J SM J SM 

Since G{u) = p — XV (u) we derive 

/ HuVG\ud^ = —k / p'^d^ + k \V{u)pd^. 

J SM J SM J SM 

Combining the last equality with yields 



{2k + I) / p^dfi-2k XV{u)pdfi 

JSM JSM 

+ I {Hufdfx- [ {K -H{\) + \^){Vufd^i = Q. 

JSM JSM 



'SM JSM 

We may rewrite this equality as: 

2 



/ (VWTIP-'-^) d, 

- I (k- HiX) + A^^^r^) iViu)r d/i + / {Huf rf/. = 0. 
JsM V 2k + l J JsM 

Combining this equality with the hypotheses we obtain Hu = 0. Note that 

K - H{X) + <K- H{X) + X^[{k + if /{2k + 1)] < 0. 
Using Hu = in (jnHI) we obtain G\{u) = p = 0. 



□ 
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Remark 4.7. If p{x,v) = qx{v,...,v) where g is a symmetric fc-tensor, then the 
condition G(p) + HV{p)/k = is just saying that q has zero divergence. For such a 
p and = 1 it suffices to assume that is Anosov [Hj. It is unknown if the Anosov 
hypothesis is enough for k > 2. The problem is open even for geodesic flows. We 
refer to j34| for partial results in this direction when k = 2. 



We begin with the following basic open problem (also raised by Wojtkowski in [37]): 

Let (f) he an Anosov Gaussian thermostat on a closed n-manifold with n > 3. Is it 
true that cf) is transitive? 

When n = 2, a result of Ghys [121 ensures that is topologically conjugate to the 
geodesic flow of a metric of constant negative curvature and thus (f) is transitive. If 
the weak stable and unstable bundles of (f) are transversal to the vertical subspace 
and M supports an Anosov geodesic flow, then a related result in TB] also shows that 
(f) is transitive. Recall that the vertical subspace V at (x, v) G SM is the kernel of 
d7T(^x,v) '■ Ti^x,v)SM — >■ TxM. Thus it is natural to ask: 

Let (f) he an Anosov Gaussian thermostat on a closed n-manifold with n > 3. Are 
the weak stahle and unstahle hundles of (p always transversal toV? 

For n = 2 this is proved in |6j but the proof requires to know apriori that is 
transitive, so both questions are intimately related. 
We now make the following useful observation. 

Proposition 5.1. Let cf) he a transitive Anosov thermostat on a closed manifold. 
Then is homologically full, that is, every homology class in ifi(5'M, Z) contains a 
closed orhit of (p. In particular is weak-mixing. 

Proof. We would like to use Theorem 1 in jSH] which gives several equivalent con- 
ditions for a transitive Anosov flow to be homologically full. They all imply that 
is weak-mixing. The one that we will use is the existence of a Gibbs state /i with 
zero asymptotic cycle As a Gibbs state we take the measure m of maximal en- 
tropy. Since an Anosov thermostat is reversible via the flip f{x, v) = {x, —v) and the 
measure of maximal entropy is unique, we see that f^m = m. 

An easy argument with the Gysin sequence of the sphere bundle vr : SM M 
shows that n* : H\M, R) H\SM, M) is an isomorphism, so given c e H^{SM, M) 
let us write c = [ti*uj] where a; is a closed 1-form in M. Since t:*uj{Ge) = oj we have: 



If the 1-form 9 dual to E is closed (but not exact) the results in j2H| assert that is 
conformally symplectic and that the weak stable and unstable bundles are Lagrangian 
subspaces. In this case one can consider the action of c/0 on the bundle of Lagrangian 
subspaces and using the Maslov cycle. Proposition 15. II and arguments similar to those 



5. Final remarks and open problems 




But f^m = m and uj o f = —uj^ hence $m(c) = for all c. 



□ 
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in [221 Chapter 2] and [6', one can show that V is transversal to the weak bundles. 
(Details of this will appear elsewhere.) Of course, in this case we already know that 
the entropy production is positive, but the transversality property may be of help in 
understanding the cohomological equation in general. 

Besides transitivity and transversality of the weak bundles with V, the other out- 
standing open problem is this: 

Let (f) he a transitive Anosov thermostat on a closed n-manifold M with n > 3. Let 
d he a smooth 1-form on M. Suppose u is a smooth solution of 

Ge{u) = ^. 

Is it true that is exact? 

For n = 2 this is proved in [H] and Theorem 14.31 provides an affirmative answer 
under certain curvature condition. 
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